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Scattering in a Simple 2-d Lattice Model* 
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Liischer has suggested a method to determine phase shifts from the finite volume dependence of the two-particle 
energy spectrum. We apply this to two models in d=2: (a) the Ising model, (b) a system of two Ising fields with 
different mass and coupled through a 3-point term, both considered in the symmetric phase. The Monte Carlo 
simulation makes use of the cluster updating and reduced variance operator techniques. For the Ising system we 
study in particular O(a^) effects in the phase shift of the 2-particle scattering process. 
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1. INTRODUCTION 

Most of the particle entries in the annual review 
of particle properties of the particle data group re- 
fer to resonances. Usually they are created in two- 
particle collisions and observed in corresponding 
cross sections and scattering phase shifts. How- 
ever, in quantum field theory unstable states still 
are one of the tantalizing problems. In the lattice 
approach, up to now efforts have concentrated on 
the identification of stable particle states. The 
determination of phase shifts and resonance pa- 
rameters has been reserved for a future of more 
powerful computers and there are only very few 
(courageous) contributions ||l|- Q. 

Recently, Liischer [gj has refined earlier ideas 
(iJI^ on how to determine phase shifts in an elas- 
tic two-particle system from the energy spectrum 
in a finite volume and together with Wolff Q the 
power of the method was demonstrated in the 
d=2 0(3) model. Here we discuss some of our re- 
sults Q where we use this method to determine 
the phase shifts in a simple d=2 model with two 
particle types and resonating phase shifts. Other 
contributions to this conference present first re- 
sults for d=3 § and d=4 § (cf. also @) where 
life is harder and statistics scarcer. The case d=2 
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provides an excellent testbed to study some of the 
juicy details of the approach. 

1.1. Method 

Let us briefly review the idea. Consider the 
scattering of two identical particles of mass m in 
a box of finite spatial extension L. The time ex- 
tension is assumed to be sufficiently large, not 
to contribute to finite size effects. The size of 
the spatial volume L and the periodic b.c, how- 
ever, are responsible for the quantization of the 
momenta. In the elastic regime 2m < W < Am 
(or 3to, depending on the theory) the allowed mo- 
mentum values k are, in the two-dimensional case, 
related to the scattering phase shift via the quan- 
tization condition 



25{kn) + fc„i = 2n7r , n e N. 



(1) 



Assuming vanishing total momentum (CMS) the 
total energy of the 2-particle state is just twice the 
energy of the back-to-back single particle states 



Wn = 2 Vm2 -f kl 



(2) 



Thus given m and a couple of low lying energy 
levels Wn one may obtain values of the (infinite 
volume) phase shift at the corresponding values 
kn ■ Varying L one may cover a whole range of mo- 
mentum values. Relation (|l|) holds in that simple 
form for d=2 but can be generalized to higher 
dimensions 0]. 



One has to take care of the foUowing restric- 
tions. 

• The interaction region and the single par- 
ticle correlation length ought to be smaller 
than the spatial volume. 

• Polarization effects due to virtual particles 
running around the torus should be under 
control. 

• Lattice artifacts will turn up in 0{a^) cor- 
rections. 

• For the determination of the energy spec- 
trum one should consider correlation func- 
tions of sufficiently many observables. 

In d=2 all these can be controlled. 

1.2. Model and Simulation 

We choose a model, where two light particles 
if couple to a heavier particle rj giving rise to res- 
onating behaviour. The action is given by 
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The values of the fields are restricted to {+1, —1}. 
The sums run over all sites (a;o,a;i) of the eu- 
clidean L x T lattice A with periodic boundary 
conditions. The 3-point term was introduced in 
a nonlocal but symmetric way, because ip^ = \. 
For g — Q this is just a system of two indepen- 
dent Ising models, each describing in the scaling 
region interacting bosonic fields with mass m{K). 
The corresponding masses have been adjusted to 
m^ ~ 0.19 and rn^ ~ 0.5 (m^ defined by the reso- 
nance peak position) by calibrating the couplings 

l^/n anCl l^'n • 

When kinematically allowed, the term propor- 
tional to g gives rise to transitions like r/ -^ 
ipip rendering rj a resonance in the ipip channel. 
We study the model at g = 0, 0.02 and 0.04. 
Throughout this work we use T — 100; the spatial 
extension L varies between 12 and 60. For each 



set of couplings and lattice size we performed typ- 
ically 2 X 10^ measurements. Our Monte Carlo 
simulation utilizes the cluster updating method 
introduced for the Ising model in iQ. The sta- 
tistical errors are estimated with the Jackknife 
method. Details of the simulation technique and 
the phase diagram can be found in l^]. 

2. OBSERVABLES 

2.1. Single particle state 

The operator of a 1^9 state with momentum 
Pi^i, = 2'Ki'/L, V = — L/2 -I- 1, ... , L/2, is given 
through 



7 X! ¥'::>;o,a:iexp(zxipi,^). 



(4) 
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where kx^ denotes a timeslice of A. Its connected 
correlation function over temporal distance t de- 
cays exponentially oc exp (— po,i/^) defining po.;^; 
in particular po,i'=o — rriip. 

For the determination of the energy spectrum a 
precise knowledge of the single particle mass and 
related finite size effects is important. We find Q, 
that our results for p^^y follow with high precision 
the spectral relation for the lattice propagator of 
a Gaussian particle with mass m, 



Pq.u = cosh (1 — cospi^jy + coshm). 



(5) 



This expression deviates from the continuum dis- 
persion relation (d.r.) po.y = {1^'^ +Pi ,y)^ by a 
leading correction 0{{ap)^). 

The observed mass, as compared to the "real" 
mass at vanishing lattice spacing and infinite vol- 
ume, is also affected by polarization due to self 
interaction around the torus. We confirm this 
behaviour and find good agreement with the ex- 
pected exponential decrease fill]. We also de- 
termined the wave function rcnormalization con- 
stants for the fields Q. 

2.2. Scattering sector 

We consider operators with total zero momen- 
tum and quantum numbers of the rj, 
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Figure 1. Phase shifts for the Ising model for two 
values of the mass, determined with the contin- 
uum d.r. (2) ; full line: S = —tt/2, dashed line: 
2-particle threshold. 






with pj — 



ipj{xi-yi) 



J-2,3...(7) 



and measure all (connected) cross-correlations 
Mnm{t) ~ {Nn{t)Nm{0))c {t deuotcs the sepa- 
ration of the time slices). The operators iVj>i 
describe two (/5-particles in the CM system with 
relative momentum 2pj . Because of the interac- 
tion they do not correspond to eigenstates of our 
model. Indeed they are eigenstates of the Gaus- 
sian model of free bosons. However, if the set 
is complete, a diagonalization of the correlation 
matrix provides the necessary information on the 
energy spectrum of this channel. The transfer 
matrix formalism yields the spectral decomposi- 



Figure 2. Like fig. 1, but now using the lattice 
d.r. (9) 



tion 
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where Vn' — {l\NnO) are the projections of the 
states |A^„0) (generated by the operators Nn out 
of the vacuum) on the energy eigenstates (^| of 
the scattering problem. 

The number of operators considered should be 
chosen larger than the number of states in the 
elastic regime 2m^ < W < 4m<^ and not larger 
than L/2 to be linearly independent. A larger set 
provides a better representation of the eigenstates 
but enhances the numerical noise. 

We work with between 4 and 6 operators de- 
pending on L and distances t = 1 ... 8. Details 
on the determination of the eigenspectrum in the 
scattering sector and on the representation of the 
eigenstates are discussed in g. 



3. RESULTS FOR THE PHASE SHIFTS 

3.1. Resonance 

In [pj we present our results for the observed en- 
ergy levels and the resulting phase shifts. Since 
the pure Ising model has an S-matrix equal to —1 
|12], the phase shift starts with — 7r/2 and then 
shows, for non- vanishing g^ a clear resonance be- 
haviour manifesting itself in a fast increase by tt. 
It can be nicely approximated by a standard ef- 
fective range resonance formula 



3.2. 0{a^) corrections 

Let us consider the Ising case {g = 0) in more 
detail. Here we increased the statistics signifi- 
cantly (1.5 million measurements) and repeated 
the calculation for m = 0.19 and a higher value 
of the mass m = 0.5 (fig.0). As mentioned we 
expect a phase shift of — 7r/2 and the results are 
in agreement with this. At small k the energy is 
close to the 2-particle threshold and the statisti- 
cal error of the energy transforms via (0) into a 
relatively larger error of k and thus 5{k). Higher 
k, on the other hand, stem from large values of 
the energy with intrinsically larger statistical fluc- 
tuations. However, due to the enhanced statistics 
we do identify a systematic deviation from — 7r/2 
increasing with k and m. We attribute this be- 
haviour to 0{a?) corrections. 

As mentioned earlier, the d.r. (g) gives the 
total energy of the asymptotic 2-particle state 
which (under the assumption of localized inter- 
action region) is just twice the energy of the out- 
going particles. Now the 0{a?) corrections of the 
single particle d.r. have been nicely described by 
replacing the continuum d.r. by the lattice rela- 
tion (^). We therefore replace (0) by the corre- 
sponding lattice expression. 



Wn — 2 cosh ^ (1 — COS kn + cosh to). 



(9) 



Our data for Wn now produce slightly different 
values of fc„ and 5 exhibited in fig.y, in better 
agreement with a constant value of — 7r/2. We 
conclude that the leading O(a^) corrections can 
be expressed by replacing the continuum d.r. by 
the lattice relation, at least in the 2 dimensional 
Ising Model. In general we suspect that lattice 
artifacts in the phase shift can be diminished by 



studying carefully the dispersion relation of single 
particle states. 

4. CONCLUSION 

We have determined phase shifts in the Ising 
model and resonating phase shifts in a model with 
two types of particles and a three-point coupling. 
We find that Liischer's suggestion for determin- 
ing these phase shifts is indeed a very reliable 
method, at least in d=2. The leading 0{a'^) ef- 
fects in the 1- and 2-particle channel may be ex- 
plained by the differences between lattice vs. con- 
tinuum dispersion relations. 
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